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Abstract 



The spontaneous decay rates of an excited atom placed near a dielectric 
cylinder are investigated. A special attention is paid to the case when the 
cylinder radius is small in comparison with radiation wavelength (nanofiber 
or photonic wire). In this case, the analytical expressions of the transition 
rates for different orientations of dipole are derived. It is shown that the 
main contribution to decay rates is due to quasistatic interaction of atom 
dipole momentum with nanofiber and the contributions of guided modes are 
exponentially small. On the contrary, in the case when the radius of fiber is 
only slightly less than radiation wavelength, the influence of guided modes can 
be substantial. The results obtained are compared with the case of dielectric 
nanospheroid and ideally conducting wire. 

42.50.Pq, 33.50.-j, 32.80.-t, 12.20.Ds 
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I. INTRODUCTION 



At first sight, spontaneous emission is a purely atomic process. However as was first 
pointed out by Purcell [1], a resonant cavity can change the decay rate substantially. At 
present it is known that not only resonant cavities but also any material body can influence 
the decay rate of spontaneous emission [2] . Moreover, the control of decay rate of sponta- 
neous emission is widely used in practice of making the new efficient sources of the light 
[3]. 

In spite of a qualitative understanding of the infiuence of material bodies on the sponta- 
neous radiation of atom, only the influence of plane or spherical interface on decay rate has 
been elaborated in detail [4-10]. 

However, the influence of dielectric flber or metallic wire on decay rate of a single atom 
is interesting from a theoretical and practical point of view. This is due to the fact that the 
charged wires are successfully used to control atom motion [11,12]. The cylindrical geometry 
is also important for investigation of fluorescence of substances in submicron capillaries 
[13,14]. A very important application of the decay rate theory in the presence of dielectric 
flber is the photonic wire lasers [15,16]. Finally, cylindrical geometry appears naturally when 
considering carbon nanotubes [17,18]. 

The influence of ideally conducting cylindrical surface on decay rates is well investigated 
both for atom inside a cylindrical cavity [19-21] and near an ideally conducting cylinder [22]. 

The spontaneous emission of an atom in presence of a dielectric, semiconductor or metal- 
hc cyhnder is a more complicated process. First investigation of the decay rate of an atom 
placed on the axis of dielectric flber was undertaken within classical approach in [23-24]. 
Recently that problem attracted new interest [25-28]. In [29-30] spontaneous emission near 
carbon nanotubes was considered. The general line of novel papers was a wide using of 
numerical methods. Unfortunately such an approach cannot answer many qualitative ques- 
tions. In this paper, we re-investigate the influence of a dielectric cylindrical surface on rates 
of dipole transitions, using analytical and asymptotic approaches. For brevity we will con- 
sider only one channel of decay. For example it may be 2P-1S channel or any other channel 
of decay. To obtain clear analytical results we will investigate flber with radius small in 
comparison with a wavelength. We will pay special attention to different mechanisms of the 
decay rate (radiative, waveguided, and nonradiative) . 

The structure of the rest of the paper is as follows. In Section II , we present general 
theory of decay rate in the presence of any body. We shall also simplify general results 
in the case of nanobodies. In Sect. Ill, we consider the decay rates for an atom placed 
in a close vicinity to a nanoflber with any complex dielectric permittivity, when one can 
use quasistatic approximations. We obtain simple analytical expression for radiative and 
nonradiative decay rates for any orientation of dipole momentum of atom. In Sect. IV, 
we consider the full clcctrodynamic problem of dipole decay rate near a nanoflber with any 
complex dielectric permittivity. We build the analytical expression for decay rate through 
contour integral in complex plane of longitudinal wave number. Then we transform general 
contour integrals to separate contributions of guided and radiating modes. In Section V, 
general expressions obtained in Section IV are applied to z and </?, p orientations of dipole 
momentum of an atom placed at the surface of lossless flber. In Section VI, wc present 
graphical illustrations and discuss the results obtained in previous Sections. Geometry of 
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the problem under investigation is shown in Fig.l. The dielectric permittivity of a cylinder 
is £, the dielectric permittivity of an outer space is equal to 1. 



II. DECAY RATE OF AN ATOM IN THE PRESENCE OF AN ARBITRARY 

BODY 

In general the spontaneous emission of an atom placed in the vicinity of any body is due 
to three different factors. First, the excitation energy can be emitted in the form of photons 
moving to space infinity. It is natural to name the corresponding part of decay rates as 
a radiative decay rate and to designate it as Y"''^^'^^^'"'^ . Second, the excitation energy can 
be transformed to photons locaUzed near or inside a body, that is, to guided modes. We 
will denote the corresponding part of decay rates as Finally, in the case of complex 

dielectric permittivity the excitation energy can be transformed into thermal energy of a 
body. We will denote a corresponding part of the decay rates as 'y'^o'^^O'd.iative ^ 

Thus the total decay rate of an atom placed near any body can be represented in the 
following form: 

^total ^radiative , ^,guided , ^nonradiative /-i \ 

In the case of excited molecules or quantum dots the total decay rate includes also the 
internal nonradiative transitions. Such transitions have no resonant nature and the influence 
of a body on these transitions seems to be insignificant. In this paper we do not take such 
transition into account. 

Let us now consider the radiative decay rate in more detail. Within classical approach 
the excited atom can be described as a linear oscillator, whose dipole momentum do is 
proportional to matrix element of dipole momentum operator. The oscillation frequency 
Uq is equal to transition frequency. The radiation power of classical dipole in free space is 
described by a well known expression [31] 



C, T-\ \ radiative 
dt stt y V ; V ; r-.oo 3 



r^dQ ^^\dof (2) 



where integration is over sohd angle, E(°),H(o) are electric and magnetic fields in free space, 
k — ^ stands for the wave vectors in free space, and c is the velocity of light in vacuum. 

If one puts any body near excited atom, the radiation power will be changed and will be 
described by the formula 



'dtj 



±. J I (e(°) + E(«)) X (h^'^) + H(^)) ll^ydn (3) 



where E^^^ and H^^^ are the electric and magnetic fields refiected by a body. 

If, additionally, the size of body (more precisely, the characteristic size of the region with 
nonzero polarization) is small in comparison with radiation wavelength and atom is placed 
near a nanobody, the radiation of a whole system will be of a dipole type. As a result, we 
will have, instead of (3), the following expression 
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radiative ,4 

'd,.,f (4) 



3 

where dtot stands for total dipole momentum of the whole system. 

As the radiative decay rate is proportional to radiation power, the relative radiative 
decay rate can be presented in the form: 



^.radiative I J 1 2 

7 I ^tot I 



70 I do 



2 



(5) 



where 70 is the decay rate of an excited atom in free space. 

Within quantum mechanical approach the radiative decay rate is described by the Fermi 
golden rule 



Y-^^-^^^^oc\doE,ac{r')\'pF{u;) (6) 

where do is the matrix element of dipole momentum operator, E^^c (r') is the strength of 
electric field of emitted photon at the atom position, and Pf ((^) is the density of the final 
photon states. Note that E^^c (r') is the solution of Maxwell's equations in free space. When 
the nanobody is present, the radiative decay rate is described by the Fermi golden rule again 



Radiative ^ (j./) |2^^ ^^^j 

but now E (r') is the solution of Maxwell's equations, which are modifies by the presence of 
nanobody. The density of state pp (cu) is still independent of nanobody. Sometimes quantity 
■l^SjM^lLpp (^^j^ is referred to as the radiative local density of states. That quantity depends 
on the presence of nanobody. 

As E„ac (r) is nearly uniform on the scale of nanobody, the nanobody acquires the dipole 
momentum 6d = aE^ac, where a is the nanobody polarizability tensor. As a result the 
influence of nanobody on solution of Maxwell's equations can be described by the following 
expression 

E = E„„c + SE, SE = GSd (8) 

where SE is the electric field of dipole with momentum Sd, and G is the tensor Green 
function of dipole source. Using symmetry of G and a, it is possible to show that 

doE{r') = dtotalKac{r') (9) 

where dtotai = do + aGdo is the total dipole momentum of classical system where matrix 
element do is a dipole source. Now substituting (9) into (7) and averaging over polarizations 
of emitted photons we again obtain expression (5) for the relative radiative decay rate. 

Thus, to find the radiative decay rate of an atom in the presence of nanobody it is 
sufficient to solve quasistatic problem and to find the total dipole momentum of the whole 
system. Let us stress once more that (5) is valid within classical and quantum electrody- 
namics. 
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In the case of waveguided photons, the Fermi golden rule remains valid, but it is impos- 
sible to obtain general expression for the case of nanobodies, because the space structure 
of waveguided photons depends substantially on the nanobody geometry. The subsequent 
analysis (Section V) shows that the decay rate into guided modes decreases exponentially 
with decreasing of nanobody size. 

Let us now consider the nonradiative losses, that is, losses due to presence of imaginary 
part in dielectric permittivity (nonzero conductivity). In this case one should use general 
expression for total decay rate of an atom near any body [33]: 

7*"*"' ^^3 doE(^)(r^a;o) 

= 1 + -Im -2— (10) 

Within this approach it is sufficient to find refiected field at the atom position E^^) (r', uo). 
It is very important that this expression is again valid within classical and quantum elec- 
trodynamics [34] - [39]. 

In general it is very difficult to find an analytical expression for reflected fleld. However, in 
the case of nanobodies one can use perturbation theory over wave vector k (long wavelength 
perturbation theory) [40]. Within such theory the electric field can be presented by a series 

doE(^y,a;o) ^ ^ ^ ^^^2 ^ .^^p ^ (11) 

"0 

where the coefficients Oi, bi, Ci, and di are determined by solving some quasistatic problems 
[40]. It is important to note that the first three terms are due to near fields, while radiation 
fields appear only starting with the fourth term, which is proportional to k^. In the case of 
a medium with losses, all of the coefficients Oi, bi, Ci, and di are complex. Now substituting 
(11) into (10) we obtain the series for the total decay rate: 

total 3 / bi Ci\ ^ 3^ , 

For nonradiative part of decay rate from (1) we have : 

^nonradiative ^total ^radiative ^guided (13) 

Due to the fact that the radiative and guided decay rates begin with the fourth term of 
expansion, the nonradiative decay rate can be presented in the form : 

^.nonradiative o / „ \ 

^Inr^ + ... (14) 



70 2 VA;3 

Thus, to find the nonradiative decay rate it is sufficient to find quasistatic refiected field. 
Let us stress once more that (14) is valid within classical and quantum approaches. 

By and large, to find radiative and nonradiative decay rates in the case of nanobodies it 
will suffice to solve quasistatic problem and to find the refiected field and the total dipole 
momentum of the whole system. Of course, the decay of an excited atom into waveguided 
modes can not be described by quasistatic approximation. 
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III. QUASISTATIC ANALYSIS OF DECAY RATE NEAR NANOFIBER 



Let us consider now a decay rate of an excited atom placed near a nanofiber. When 
the distance between the atom and fiber, and the radius of fiber are substantially less than 
radiation wavelength, the fiber is polarized only near the atom. 

As a result, the radiation of atom + fiber system will be of the dipole-type and can be 

described by (5). 

The total dipole momentum dtot can be found from the solution of quasistatic problem 

7T °. (15) 

where the charge density pc at the point r is derived from the dipole momentum of atom 

pc^-{doV)d^'\r-r')e-'^\ (16) 

S^^^ (r) is the three-dimensional Dirac delta function and V' means gradient over radius vector 
of the atom, r'. Hereafter we will omit the time dependence of the fields. The continuity 
conditions for the tangential components of E and the normal components of D at the 
surface of a cylinder with dielectric constant e should be provided as well. 
Introducing a potential ip by 

E = -V(doV')<^(r,r'), (17) 

we obtain, instead of (15), the Poisson equation. 



—V^ip = AttS^^^t — r'), outside cylinder 
— V^<^ = 0, inside cylinder 

It is convenient to represent the solution of problem (18) in the form 

if = ipQ + ip2, outside cylinder 
(p — (fi, inside cylinder 

where (fo is the free-space potential given by 



(18) 



(19) 



'fo = (20) 
|r — r I 

This free space Green function can be expanded in cylindrical coordinates system (p, ip, z) 
(Fig.l) in the following series [41] : 



1 _ 2 
— r'l TT 



oo °? 

^ (2 - 5m,o) cosm ((/? - if') / dh COS h{z - z')Km {hp') Im (hp) {p < p') (21) 

m=0 n 



where 5mfi is the Kronecker symbol, K and I are the modified Bessel functions [42]. 

Using this expansion we can found expressions for (pi and <^2 by usual mode matching. 
For potential outside fiber we have the following expression (p > p') 
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2 oo °? 

= - E (2 - C) cosmic - / dh cos h{z - z')Km {hp') Km {hp) Gm {ha) (22) 



2 °° 

m=0 

where reflection coeflicients Gm look as follows: 



Gm{s) 



(s-1) 



dlm(s) 



ds 



{s) 



^I^{s)-s'-^Km{s) 



(23) 



It is important to note that for negative values of dielectric constant, Re{e) < — 1 
(metals), the denominator of (23) is equal to zero for some values of integration variable h. 
It means that under such conditions, the guided modes do occur in the system. 

To determine the total dipolc momentum of the system one should find the asymptotic 
value of (22) at large distances, + 2;^ 00. In the dielectric case (Re(£) > 1) the main 
contribution to (22) is due to m=l term and has the form 



¥>2 



-—cos {ip 
1 P 



P 



Comparing the electric potential of reflected fleld 

(^f = (doV) ^2 

with dipole potential 

~ _ dtotr 

Vdip — 5~ 



(24) 



(25) 



(26) 



one can flnd the dipole momentum of flber, 5d. In the case of p-orientation of dipole 
momentum we have 



Sdp = do,, 



la' 



while in the case of (^-orientation the dipole momentum of flber will be 



Sdip = -do,( 



e-la' 



e+lp' 



,/2 



(27) 



(28) 



In the case of z-orientation the dipole momentum induced in flber is equal to zero. 
Now combining the dipole momenta of nanoflber and atom and substituting the result 
in (5) we obtain the radiative decay rates for an atom near a nanoflber 
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radiative ' 



7o , 

radiative ' 



1 + 



7 



70 , 

radiative ' 



1 - 



£ + 1 p'2 

- 1 a2 



e+lp 



/2 



(29) 



7 



= 1 



7o 
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In the case of an atom placed at the fiber surface these results agree with those obtained 
in the case of a prolate dielectric nanospheroid [43,44]. Besides, these results agree with 
decay rates of an atom placed on the axis of dielectric fiber [23,24]. Indeed, for p' = a we 
have from (29) 



1 



radiative " 



7o 



2e 



e + 1 



7 



radiative ^ 



7o 



e + 1 



while for the atom at the nanofiber axis one has [23,24]. 
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radiative ^ 



7o 



(30) 



7 



radiative ' 



7 



radiative ' 



7o 



7o 



e + l 



7 



radiative ' 



= 1 



7o 



(31) 



The only difference is for p- oriented dipole, where we have enhancement by e^. The 
physical interpretation is as follows. From QED point of view, decay rate is due to inter- 
action of the dipole with electromagnetic modes, modified by a fiber. The tangential {if, z) 
components of mode electric field are continuous. So the decay rates of dipoles with (99, z) 
orientations at the surface and in the body of the fiber should be the same. On the con- 
trary, the p - component of mode electric field is discontinuous at the surface (the normal 
component of = eE is continuous). This explains the difference between decay rates of 
p-oriented dipole inside and outside nanofiber. 

We should also mention that our results do not agree with decay rates of an atom inside 
photonic wire found in [15,16]. We believe that results [15,16] are misleading because of a 
very crude approximation of fiber by plane waveguide. 

In the case of nanowires, that is, in the case of Re{e) < —1, the excitation of guided 
modes is possible and one should add additional (very important !) terms to (29). The 
detailed investigation of influence of guided modes on decay rates near metallic nanowire 
will be presented elsewhere [45]. 

Analogous additional terms occur in the case of dielectric fibers and they are also due 
to the guided modes. In following Sections we will show that those terms are exponentially 
small for nanofibers. As a result one can use quasistatic formulae for the nanofiber safely. 

The results obtained, (29), are vahd only for dielectric and metallic nanocylinders. The 
situation is dramatically changed in the case of an ideal conductor (|£| — > 00) , when reflection 
coefficients become equal to 



G. 



ideal conductor 



is)- 



Imjs) 



(32) 



and main contribution to far field is due to m = term. As a result the potential of refiected 
field for p-orientation of dipole decreases at infinity more slowly than (26). It means that 
dipole momentum induced in an ideally conducting nanowire tends to infinity when wire 
radius goes to zero. Respectively, the decay rate tends to infinity too. This fact is in 
agreement with the exact solution of the elcctrodynamic problem [22], where it was shown 
that the decay rate of radially oriented dipole tends to infinity when cylinder radius tends 
to zero. 

Up to now we considered the radiative decay rates only. However even for dielectric 
nanobodies the nonradiative losses can be substantial. Here the nonradiative losses are due 
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to losses inside nanobody, which are proportional to imaginary part of dielectric permittivity 
(or conductivity). To take the nonradiative losses into account, one should use (14). 
In our case, the quasistatic solution E^^^ is given by 

E(«) = -V (doV) (^(2) (r, r') , (33) 

where <^*^^) (r, r') is given by (22). Substituting this expression into (14) we obtain the 
expressions for the nonradiative decay rates near a nanofiber: 

/ ^nonradiative \ o oo °? 

-n=0 

^ — / dhKl {hp')Im (Gm (ha)) (34) 



7o 



' ^nonradiative \ g oo ^2 

70 tt/c^JoP' 



' ^nonradiative ^ 



7o / ^ T^k^ 



oo ^ 

5^ (2 - 5^,o) / dhh^Kl (hp') Im {Gm (ha)) 

m=0 i 



If an atom is placed in close vicinity to surface of a nanofiber {p' a), the above ex- 
pressions can be reduced to 



' ^nonradiative \ 

— Im 



7o Jp Ve + iy 16A;3(p'-a) 



nonradiative \ / ^.nonradiative \ / tr 1 \ Q 

(35) 



' ^nonradiative \ I ^nonraaiative \ / _ i 

Im 



V 70 V 70 J, \e + IJ 32k^p' - ay 

It should be noted, these expressions are similar to the case of an atom near plane 
dielectric interface, where the reflected field can be described by image dipoles. Comparing 
(35) and (29) one can see that for usual dielectrics with low losses (for fused silica Ime ~ 
10"^) the nonradiative losses are very small for any reasonable distance of an atom from a 
surface. However, when the atom is placed very close to the surface, the nonradiative losses 
can be enhanced substantially. 



IV. FULL ELECTRODYNAMIC APPROACH TO DECAY RATES NEAR 

DIELECTRIC FIBER. 

In previous Section the decay rates were found within quasistatic approximation without 
taking into account guided modes. In the present Section, within full Maxwell's propagation 
theory, we will find the exact expressions for decay rates, which include the contributions 
from guided modes. 

According to [34] - [39] , classical and quantum-electrodynamic calculations give the same 
results for dipole transition rate normalized to its vacuum value. In present section we shall 
investigate the influence of a dielectric cylinder on transition rates within classical approach, 
where the full decay rate can be expressed through classical reflected fleld at the atom 
position (14). Thus the reflected field must be found to determine the total decay rate. To 
find the reflected fleld, it is necessary to solve the full system of Maxwell's equations where 
dipole momentum of oscillator is a source and to use appropriate boundary conditions. 



9 



To find the reflected field we follow an approach develoded by Katsenelenbaum [23,24] and 
Wait [46]. According to that approach one should expand all fields over cylinder harmonics. 
For longitudinal components we have the following expressions: 

m=—oo m=—oo 

oo 



m=— oo 
oo 



m=—oo m=—oo 



(36) 



m=— oo 

oo 



m=— oo 

BP= E Jdhe'-^'^+'^'BfXhip)- § Jdhe'^'P+'^'JmMdmh 



Here superscripts R and T and subscript correspond to dipole field in free space, 
reflected and transmitted inside cylinder flelds, respectively, and vi = \J ek^ — h? and 
V2 = y/k'^ - h? are the radial wavenumbers inside and outside fiber. The 

(^mhi bmhi C-mhi dmh 

coefficients are to be determined. We choose the branch cut such as Imiyi) ,Ifn{^2) > 
in complex plane of longitudinal wavenumber h. To ensure the decreasing of fields at space 
infinity (p — > oo) the integration in (36) should be over path Ci, which is shown in Fig. 2a. 

The rest of the field components (p,^?) can be expressed through z-components of electric 
and magnetic fields: 

P ih_ dEz^mh km 



(37) 



rp hm rp ik_ dBz „^i^ 

^tp,mh — pi,2^z,mh ^2 Qp 

TD kem Tp 4_ ih. ^^'^t^h 

J->p,mh — pj,2 ^z,mh "r ^2 Qp 

TD ike d^z,mh mh td 

J->tp,mh — 1,2 Qp pj^2J->z,mh, 

In (36), (37) the subscripts m and h denote the appropriate Fourier transformation over 
(f and z, and v — Vi or v-i for corresponding space region. 
For free field one can obtain the following expressions [31]: 

B = rotA, E = \rom 

gifelr-r'r (38) 

A = —ikdiQ- - 

|r — r I 

where r = [p^ip^z) and r' = (p', (^', z') are radius vectors of observation point and atom 
position. 

Using the well known expression [47] (p < p') 



Y: f dh e'^^^-^'^+'^^^-^'^mM H^^ (u,p') (39) 

m=— oo^^ 
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where integration part should be as shown in Fig. 2a, one can find the expressions for longi 
tudinal components of the free dipole fields near the fiber surface. For p oriented dipole ii 
the region between fiber and dipole {a < p < p') we have 



B,.,mH = -'-^Jm M Hg) M 



2p' 

Z=V2p' 



(40) 



B,,,raH = ^ {^,Jm M Hg) ^ + fj^ {z) ^^^^^j-Hg) (.)_^^,) 

while for z oriented dipole we have 

Boz,'mh = 

Eo.,mh = '^Jm M Hg) {U^p') 

Eo^,mH = -'^Jm M Hg) M (41) 

Finally, for (f orientation of dipole momentum we have the following expressions 

Bo„mH = ^ {Jm M i-).^.., + fHg^ M (.)|_^J 

E,.,mH = --^Jm M Hg) {y,p') 

E,„^n = ^ (k^ijm iz),^.,,£Hg' i-).^.., + jM^)^- (^^^) ^^'^ 

To find coefficients in transmitted and reflected flelds one should take into account bound- 
ary conditions on the surface of dielectric cylinder. As a result we have a system of 4 
equations for 4 unknown coefficients. 

Eg^ (^2) 0.mh - Jm (^1) Cmfe = -Eoz,mh {P = «) 
Hg^ {Z2) bmh - Jm {Zl) dmh = -Hoz,mh {p = «) 

f^Hg) (z,) a^n + f,^,Hg) {z,) b^H - ^Jm (^1) c^h - f^i^Jm (^1) d^H = Eo^,mk {p = a) 

%^Hg) {Z2) amH - fHg) {Z2) bmH - '^^Jm (^l) C^. + fjm (^l) d^h = -H,^,mh {p = «) 

(43) 



(42) 
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where we use abbreviation 2:1^2 = ^1,20. 

The reflected electric fields are determined only by a^h and bmh coefiicients, which can 
be simplified to 

_ na _ nb , , 

^mh - p2 + ' ^rnh " p2 + gp 1^4 j 

where 

na = vlvlaJm {Zl) PEo^^rnh + T^l {Jm izi)hmP+kaeUi-£^Jm (^l) Q) Eoz,mh 
^ivlvlaJm {Z\) QBo^^rnh " imhl/il/2Jm {Zl) SBoz,mh 



(45) 



nb = vlvlaJm {Zl) PBo^^rnh + T^l {jm {zi)hmP-akUi-£^Jm (zi) i?) Boz,mh 



(46) 



and 



P^hmk''Jm{z,)H(^') {Z2) (s-1) 

R = u^u^ak {u^J^ {z^) ^^H^ (^,) _ :y,eHW (z,) A.J^ [z,)) (47) 



By substituting the expressions for reflected fleld (36)- (37) into general expression (10) 
we obtain final expressions for total decay rates: 

OO r -1 

h] =i + ^ir^!!!zZ^ ll ^ (48) 

V 70 2 do^pk^ 



OO r 

E §dh^[ -ikv^j-Hg) {z) bmh - ^Hg) {z) amh 



Z=l/2p 



V 70 2 do,^k^ 

OO 

/ ^.total \ 

V—] + (50) 

As was mentioned above, the important feature of dielectric fiber is the presence of 
guided modes, which differs substantially from free space spherical (or cylindrical) waves. 
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Such modes are effectively used in optical communications lines. Note that there are no such 
modes in the case of dielectric sphere or ideally conducting cylinder. So-called whispering 
gallery modes, which occur in the case of sphere or cylinder, are the decaying ones even in 
the case of lossless materials. 

Thus there are two types of modes in presence of dielectric fiber: free space radiation 
modes and waveguided modes. From mathematical point of view these modes correspond 
to different types of a spectrum. The guided modes correspond to discrete spectrum while 
radiating modes correspond to continuous part of the spectrum. In this connection the 
expressions (48)- (50) are not fully suitable for further analysis, because the guided modes 
are not separated here. 

The appearance of guided modes is connected with poles in subintegral expressions in 
(48)- (50). One can show that in the case of a lossless dielectric these poles are situated on 
the real axis of h between k and ky/s and the number of poles is finite [23,24]. 

Now transforming integration contour from Ci to C2, as shown in Fig. 2b, and applying 
the residue theorem, one can separate radiating and guided modes 



(ytotal 
70 



1 + f Im- 



-00C2 



-|-37rRe J2 J2 Res 

m=—oo hem 



1 



(51) 



h=ha 



\ 10 J in 2 



-|-37rRe E E Res 

m=—oo ha.m 



1 



(52) 



h — ha,T\ 



(ytotal 
70 



E f'"' 

1 + llm^^'- 



-|-37rRe E E Res 

m= — 00 ha,m 



do,zk^ 



do,zk^ 



(53) 



h=ha,', 



where the sum is over all poles ho,,m of subintegral functions , Res means residue, and a^^, hmh 
are described by (44)- (47) . 

In expressions (51)-(53) the first term corresponds to spherical waves running to infinity 
and nonradiative losses, while the second term corresponds to guided modes. It should be 
emphasized that in the case of lossless media the vertical part of branch cut (along imaginary 
axis of h) gives no contribution to decay rates. For lossy media this part is very important 
because it contributes to nonradiative losses. The nonradiative decay rates found in previous 
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section (expressions (34)) are the asymptotics of integral over the vertical part of branch 
cut when ka moves to zero. 

As was mentioned above, in the dielectric case there is only a finite number of guided 
modes with longitudinal wavevectors between k and k^/e. Moreover, for small enough fiber 
(for nanofiber!) with 

ka < (54) 

(where jo,i ~ 2.4048 is the first root of Jo {z)) the only one guided mode with m — ±1 
exists. Sometimes such modes are referred to as main or principal modes. 

The dependence of longitudinal wavenumber h on cylinder radius or frequency is deter- 
mined by dispersion equation 

D = + QR = 0, m = ±l (55) 

where P , Q and R are defined by (47). 

In the case of nanofibers, ka << 1, the asymptotic solution of (55) can be presented in 
the form: 

/ 2 e+1 e+1 \ , , 

where 7 = 0.5776 is the Euler constant. 

The exact and asymptotic solutions of (55) are shown in Fig. 3 , where one can see that 
asymptotics (56) presents solution of (55) correctly if ka < 0.8. In what follows we restrict 
ourselves to the case of nanofiber, where condition (54) holds true. 

To calculate decay rates into guided modes one should know residues of the corresponding 
expression. The residues can be found if one knows the asymptotic behavior of resonant 
denominator near pole. In the case of nanofiber the denominator can be approximated by 

D^P^ + QR^ A;^ ^^ ~P {kaf {{h/k) / {h/k)^^ - 1) (57) 




where {j:) is given by (56). 

Let us stress once more that the expressions (51) - (53), (56) and (57) are valid in the 
case of any complex dielectric permittivity. In the case of metallic cylinder, that is, in the 
case where Re{e) < —1, the poles of subintegral function are also near the real axis of h. 
But now they correspond to symmetric (m = 0) guided modes. More detailed investigation 
of infiuence of symmetric (m = 0) guided modes on decay rate will be presented in a separate 
publication [45]. 



V. DECAY RATES NEAR SURFACE OF NANOFIBER WITHOUT LOSSES. 

The expressions (51) - (53) fully describe the problem of spontaneous emission of an atom 
placed near a cylinder made of any material. However, those expressions are too complicated 
to understand real picture of decay rate. Our goal is to find simple analytical expressions 
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allowing one to estimate decay rate near dielectric {e > 1) cylinder with radius, which is 
substantially smaller in comparison with wavelength, ka << 1. Moreover, to obtain simple 
asymptotes we restrict ourselves to the case of an atom placed at the surface (p' — > a) of a 
lossless dielectric cylinder. 



A. z-orientation of dipole 

Substituting the expression for exciting external fields (41) into general expression for 
decay rate (53) one can represent the decay rate for z-oriented dipole momentum in the 
form: 



fe3 



--Im E I ''^"'"^ H^\''2pTJm(zi)Q 



m=— OO 



(58) 



-3Re E E Res 

m=— oo hem 



H^\u2p'fJm{zi)Q 



where Q and D are defined by (47) and (55), and 2:1,2 = ^^1,20- 

Note, the first line of (58) coincides with decay rate near ideally conducting cylinder [22]. 
This expression is valid for real dielectric permittivities. 

In the most interesting case of an atom near surface of nanofiber (p' = a) using identity 
E Jm (^) — 1 one can simplify (58) to a more compact form 



/^totoix ^ _3j^ y 'r dh,.^4 Hi^Hz2)Jm(zi)Q 
\ -TO Jz ^ m=-ooi ^ 



-3Re E E Res 

m=-0O ha,m 



fe2 D 



(59) 



The asymptote of (59) for ka <^ 1 has the following form 



(^).-l-^(60(7 + lnH-47 



45£-15 



{kaf + 0[{kay 



+ 7i^exp 



2 £+1 , £+1 



(ka)" ^^^t- y {kafe-1 



-27 



(60) 



Here the first line describes the radiative losses while exponentially small second line 
describes contribution of the principal guided mode with m = ±1. 



'B. (f - oriented dipole 

Let us now consider the case of a dipole having (p - orientation of dipole momentum and 
being located in close vicinity to the surface of nanofiber (p' — > a). In the case of lossless 
media the general expression (52) can be simplified to 
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(^) =l-iRe § ri[^J{z,)H{z,) + H'{z,)J'{z,] 

\ / in *ri — — n 2 



m=—oo 

+ -Im y ,, h'^m'^V'2J{zi)H{z2){vlJ{zi)H'{z2)-iylJ\zi)H{z2))+u^^^ 



m=— cx) 



(61) 



+3Re E E Res 

m=-O0 ha.m 



h'^m'^V2J{zi)H{z2){vlJ{zi)H\z2)-vlJ' {zi)H{z2))+uivlkaH' {Z2)J' {zi)R 

kz2D 



where R and D are defined by (47) and (55). Here for brevity we omit indices (m) in Bessel 
and Hankel functions of first kind and use prime to denote derivative of Hankel and Bessel 
functions. 

Using the identities 



E Jl {z) = 1; E n'Jl (z) 



^2 n=OD 

E 



r / dj„ (z) \ 1 



2' n^ooK dz 



(62) 



one can show that first line in (61) is equal to zero. Finally the expression for radiative 
decay rate of a (/^-oriented dipole placed at the surface of nanofiber acquires the following 
form 



\ ^0 n ^^±1^1 kz2D 



+3Re E E Res 

m=-00 ha.m 



h?m?V2J{zi)H{z2){vlJ{zi)H' {z2)-ylJ' {zi)H{z2))+VT,ylkaH' {Z2)J' {zi)R 

kz2D 



(63) 



The asymptote of (63) for ka « 1 has the following form 

/^totai\ _ / 2 \^ I (£-l)(75£2+2081-1680(7+ln(A;a)))(fea)^ f f 1 \^\ 

l^J^-li+iJ + 300^ + UI^{ka) ) 



J 48 



(e-iyika)" 



(64) 



Here the first hne describes the radiative losses while exponentially small second hne 
describes contribution of the main guided mode m — ±1. 



C. p - oriented dipole 

The case of radially oriented dipole momentum is more complicated for analysis. So we 
again restrict ourselves to the most interesting case, when atom is near to dielectric surface, 
that is, we consider p' — a case. In the case of lossless media the general expression (51) 
can be simplified to 
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m=—oo 



2j!^? J(zi) j 2h^k-'n,-Z2(:-l)J(z0IUz2)-Il\:.2)+4h-QiI'-(^2)-i>iH^U{z2fR 



-3i?e E E Res 

m=-C!0 



Pif //(Z2) I 



(65) 



where Q,R and D arc defined by (47) and (55). 

In the case e — > oo, from integral over cut (first two hnes in (65)) one can reveal the 
decay rate for an atom near an ideally conducting cylinder [22]: 



/ total \ 6 °° r 

V 70 J, T^^mt^ooi 



(i/2a)' 



1 6 ^ } 
2 + ^ /J / , 

171.= — DC) 

I k (1/20) 



(66) 



which confirms a rather complicated algebra. 

The analysis shows that for the case of cylinder of small radius ka\J e — 1 < jo,i ~ 2.4048 
only residue from term with m = ±1 will give contribution. As a result, the asymptote of 
decay rate of atom placed at the surface of small dielectric fiber takes the following form 
ka 0: 



( ^totai\ / \2 £2(£_]^)Cj^5£2_,_gQg_,_2201-1680(7+ln(fca)))(fca)^ „ /,, ,4\ 

(V) =(iTT) + ^ + 0((H) 



p 



(67) 



Here the first line describes the radiative losses while exponentially small second line 
describes contribution of the main guided mode m = ±1. 



VI. DISCUSSIONS AND GRAPHIC ILLUSTRATIONS 

We should note before all that asymptotes (60) ,(64) and (67) do agree with quasistatic 
results (29), (30) in the limit ka 0. This proves quasistatic calculations and confirms 
complicated algebra of Sections 1V,V. 

Thus, to estimate radiative and nonradiative decay rates near dielectric nanofiber with 
arbitrary complex permittivity e (Re (s) > 1) one can use expressions (29,34). To estimate 
contribution of guided modes in nanofiber one should use generalization of expressions ob- 
tained in Section V: 

/^gmded\ / 43^2 / 2 £+1 £ + 1 \\ 
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However, the asymptotics found do not allow us to determine the region of their applica- 
biUty. To find region of applicability of our results wc calculated the decay rates according to 
full formulae and compare them with asymptotes obtained in previous Section. The results 
of comparison are presented in Figs. 4-12. 

First of all from Figs. 4,7,10 one can see that our asymptotic expansions are good 
enough for ka < OA{e — 3) . When the dielectric permittivity is increased the region of 
applicability of longwave asymptotes is reduced. For e = 10 (Figs. 5,8,11) applicability 
region is ka < 0.1, while for e = 30 (Figs. 6,9,12) applicability region is ka < 0.05. One can 
suppose that generally our asymptotics are good for 

ka<l/e (69) 

The restriction is more rigid than it may appear from cursory examination {y/eka < 1). 

This region corresponds to rather thin nanofiber. For example, for £ = 3 fiber radius 
should be about A/20! The case of large dielectric permittivity should be treated carefully 
because the limits ka —> and £ ^ oo do not commute. To investigate the case of ideal 
conducting nanowire (|£| — > oo), one should take limit of general expressions (59), (63), (65) 
and only then investigate asymptotics for A;a — > 0. The analogous situation takes place in 
the case of planar interface [48] or for atom near prolate nanospheroid [43,44]. 

In the region 1/e < ka < our asymptotics give satisfactory approximation. The 

rest of nanofiber radii. 



1/Vs <ka< 2.4/Vs- 1 (70) 

should be analyzed numerically (see Figs. 4-12). From these figures one can see that for large 
enough nanofibers (70), the decay rates increase with increasing of dielectric permittivity. 
The most substantial enhancement is observed for p and z orientations, where enhancement 
of total decay rates can reach value about 35, and 23 for £ = 30 (Figs. 6, 9). The most 
important feature of p- and z- orientations of dipole momentum is very efficient excitation 
of guided modes. On the contrary, the influence of guided modes on spontaneous emission 
of (/^-oriented dipole is rather small. 

To trace the influence of guided mode on total decay rate we plot the ratio of decay rate 
into guided modes to total decay rate 



guided 



R = SJOI (71) 

f "y \ g'^'i'died / 7 \ ''"''■'^iotive V / 

\lo) VtoJ 

In (71) we neglect a contribution from the nonradiative processes. 

This quantity is very important for determining the laser threshold [49]. From Figs. 13- 
15 it is seen that the spontaneous emission coupling efficiency (5 is rather high for p - and z 
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- orientations even in the case of monomode nanofibers! It is interesting that the asymptote 
describes the couphng efficiency well in rather wide region (Fig. 13). The dipoles with tp - 
orientations of momentum have a small spontaneous emission coupling efficiency (Fig. 14) 

Finally, in Fig. 16 the total decay rates near a nanofiber are compared with decay rates 
near an ideally conducting cylinder. From the figure one can see that for large enough 
dielectric permittivity the decay rate near nanofiber can be greater than decay rate near 
ideally conducting cylinder. Again, that effect is due to excitation of principal mode in a 
nanofiber. 

VII. CONCLUSION 

In the present paper the decay rates of an excited atom placed near a dielectric fiber 
are considered. The main attention was paid to the case of cylinder with radius which is 
small in comparison with radiation wavelength (nanofiber), ka < 2.4/ — 1. The decay 
rates arc found within quasistatic as well as full electrodynamic approaches. It is proved 
that quasistatic approximation works well for a nanofiber with ka < 1/e. In contrast to 
quasistatic solution the exact one has additional terms from guided modes, which exist 
even for nanofiber of arbitrarily small radius. However the contributions from such modes 
decreases exponentially when cylinder radius tends to zero. For large enough nanofiber, 
1/e < ka < lAj \Je — 1, the influence of guided modes on the decay rate is substantial. 

The results obtained can be useful as for estimation of decay rates and for understanding 
of interplay between different decay channels. The results obtained are in agreement with 
those for an atom placed near dielectric or metaUic nanospheroid [43,44]. 

In the present paper we pay attention to the case of dielectric nanofiber with positive 
dielectric permittivity. However, our results can be applied to investigation of decay rates 
near metallic nanowire with negative dielectric constant. In the case of nanowires the qua- 
sistatic expressions (29) remain valid for description of radiative losses, but one should add 
to them a contribution, which is due to excitation of symmetric guided modes. The detailed 
analysis of decay rates of an atom placed near nanowires will be presented elsewhere [45] . 
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FIGURES 




FIG. 1. Geometry of the problem. 
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FIG. 2. Integration path for contour integral over longitudinal wavenumber h. 
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FIG. 4. Relative decay rates of an atom with p orientation as a function of radius of fiber ka 
[Eq.(65), solid lines]) and its asymptotic expression [Eq.(67), dashed lines], (atom on the surface, 
e = 3) 
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FIG. 5. Relative decay rates of an atom with p orientation as a function of radius of fiber ka 
[Eq.(65), solid lines]) and its asymptotic expression [Eq.(67), dashed lines], (atom on the surface, 

£=10) 
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FIG. 6. Relative decay rates of an atom with p orientation as a function of radius of fiber ka 
[Eq.(65), solid lines]) and its asymptotic expression [Eq.(67), dashed lines], (atom on the surface, 
£ = 30) 
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FIG. 7. Relative decay rates of an atom with z orientation as a function of radius of fiber ka 
[Eq.(59), solid lines]) and its asymptotic expression [Eq.(60), dashed lines], (atom on the surface, 
£ = 3) 
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FIG. 8. Relative decay rates of an atom with z orientation as a function of radius of fiber ka 
[Eq.(59), solid lines]) and its asymptotic expression [Eq.(60), dashed lines], (atom on the surface, 
e = 10) 
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FIG. 9. Relative decay rates of an atom with z orientation as a function of radius of fiber ka 
[Eq.(59), solid lines]) and its asymptotic expression [Eq.(60), dashed lines], (atom on the surface, 
£ = 30) 
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FIG. 10. Relative decay rates of an atom with (p orientation as a function of radius of fiber ka 
[Eq.(63), solid lines]) and its asymptotic expression [Eq.(64), dashed lines], (atom on the surface, 
£ = 3) 
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FIG. 11. Relative decay rates of an atom with if orientation as a function of radius of fiber ka 
[Eq.(63), solid lines]) and its asymptotic expression [Eq.(64), dashed lines], (atom on the surface, 
£ = 10) 
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FIG. 12. Relative decay rates of an atom with (p orientation as a function of radius of fiber ka 
[Eq.(63), solid lines]) and its asymptotic expression [Eq.(64), dashed lines], (atom on the surface, 
e = 30). On the inset: 1 - decay rate into guided modes, 2 - asymptotic expression of decay rate 
into guided modes, 3 - asymptotic expression of total decay rate, 4 - asymptotic expression of 
radiative decay rate 
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FIG. 13. The spontaneous emission couphng efficiency [Eq.(71)] as a function of radius of fiber 
ka (solid lines ) and it asymptotic expressions (dashed lines) (p - orientation, e = 10) 
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FIG. 14. The spontaneous emission coupling efficiency [Eq.(71)] as a function of radius of fiber 
ka (solid lines ) and it asymptotic expressions (dashed lines) ( ip - orientation, £ = 10 ). 
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FIG. 15. The spontaneous emission coupling efficiency [Eq.(71)] as a function of radius of fiber 
ka (solid lines ) and it asymptotic expressions (dashed lines) (z - orientation, e = 10 ). 
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FIG. 16. Relative linewidth [Eq.(65)] of dipole transitions with p orientation as a function of 
fiber radius ka for different £ = 3, 10, 30 and for ideally conducting cylinder (atom on the surface). 
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